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SECTION  1 


INTRODUCTION 


Entrained  debris  constitutes  a  widely  recog¬ 
nized  hazard  associated  with  the  firing  of  rockets  and 
recoilless  rifles.  Igniter  wires,  nozzle  plugs,  propel¬ 
lent  chunks,  ground  debris,  or  anything  else  in  the  path 
of  escaping  propellent  gases  can  be  entrained  and  accel¬ 
erated  to  speeds  sufficient  to  injure  personel  or  equip¬ 
ment.  Characterisation  of  this  debris  hazard  may  thus 
be  an  important  facet  in  the  design  of  effective  weapons. 

Section  2  of  this  paper  presents  a  model  for 
determining  a  region  of  hazard  In  the  vicinity  of  small 
rockets  and  tube  launched  weapons.  The  model  is  developed 
in  three  parts.  The  first  considers  motion  of  debris  in 
a  flow  field  generated  by  expanding  propellent  gases. 

The  second  considers  motion  of  debris  through  ambient 
air.  The  third  considers  impact  of  the  debris  -with 
personnel  or  equipment.  These  three  parts  are  then  com¬ 
bined  to  give  relations  for  estimating  the  extent  of 
the  debris  hazard  area. 

Section  3  discusses  use  of  the  model  and 
analyzes  the  Viper  light  antitank  rocket,  described 
in  Appendix  C,  as  an  illustrative  example.  This  section 
can  be  used  Independently  of  Section  2.  The  model  is 


presented  in  this  section  both  in  the  form  of  tabulated 
equations  and  in  the  form  of  dimensionless  plots.  Para¬ 
metric  sensitivity  of  the  model  is  also  discussed  here; 
debris  drag  coefficient  is  identified  as  a  critical 
parameter.  Appendix  D  gives  computer  and  calculator 
codes  for  implementing  this  model. 

Section  ^  discusses  shortcomings  of  the  model. 
While  the  debris  hazard  area  is  estimated  with  "worst 
case"  considerations,  there  are  circumstances  in  which 
particles  can  be  projected  outside  it  with  dangerous 
velocities.  Those  circumstances  can  arise  when  particle 
shape,  plume  characteristics,  or  bacicblast  area  give 
rise  to  anomolous  particle  deflections. 


. . ip . . . . 
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SECTION  2 


MODEL  DEVELOPMENT 


Tills  section  develops  the  basic  debris  hazard  model 
for  small  rocket  backblast.  The  equation  of  motion  for  a  par¬ 
ticle  in  a  flow  is  discussed  and  shown  to  be  too  complex  for 
exact  impiimentation .  Several  simplifying  assumptions  are 
then  made  which  lead  to  upper  bounds  on  minimum  safe  standoff 
distance  and  maximum  dispersion  angle.  Ideas  from  fluid 
mechanics,  gas  dynamics,  and  a  semiempirical  skin  penetra¬ 
tion  criterion  are  incorporated  into  the  debris  hazard  model. 

,  An  effort  has  been  made  to  express  the  model  in  terms  of  readily 
attainable  parameters  and  to  keen  the  number  of  those  para¬ 
meters  at  a  reasonable  minimum. 

A  particle  entrained  in  a  gas  flow  is  acted  on  by 
aerodynamic  forces  and  by  body  forces.  Equating  those  forces 
to  the  rate  of  momentum  change  of  the  particle  leads  to  an 
equation  of  motion  for  the  particle.  The  aerodynamic  forces 
are  conventionally  resolved  into  lift  and  drag  forces  which 
are,  respectively,  normal  and  antiparallel  to  the  particle's 
motion  relative  to  the  gas.  Gravity  is  the  only  body  force 
•  of  significance  to  the  debris  hazard  problem* 

For  the  present,  we  will  ignore  lift  and  gravity. 

She  effects  of  those  forces  will  be  discussed  in  Section  4* 

Even  so,  a  particle's  equation  of  motion  is  complex* 


Y 


m 


= 

dfc 


'2 


CdAp 


-  v  !  ( v 
P'  g 


v  ) 
p 


(2.1) 


where  m  is  particle  mass,  v^(t)  is  particle  vector  velocity 
as  a  function  of  time,  Cd  is  particle  drag  coefficient,  p 
is  gas  density,  A  is  particle  velocity-wise  projected  area, 
and  Vg  is  gas  vector  velocity.  In  general,  will  be  a  func¬ 
tion  of  velocity,  gas  temperature,  and  particle  dynamics  while 
A  will  be  a  function  of  particle  rientation.  In  addition, 

p  and  v  will  have  time  varying  spatial  distributions  due  to 

s 

time  variation  in  the  gas  source,  to  turbulent  mixing  with 
ambient  air,  and  to  dynamic  coupling  with  the  particle  motion. 

The  problem  in  this  form  is  intractable.  Never¬ 
theless,  it  is  possible  to  develop  equation  (2.1)  to  arrive  at 
an  estimate  for  maximum  particle  velocity  and  maximum  parti¬ 
cle  dispersion.  Such  considerations  v/ill  suffice  for  definl-^ 
tion  of  an  area  of  debris  hazard. 

Four  assumptions  will  greatly  simplify  analysis 
while  retaining  the  essential  features  of  the  debris  hazard 
problem.  First,  a  particle  will  reach  a  higher  velocity  in  a 
fully  developed  plume  than  in  a  partially  developed  one. 
Second,  a  particle’s  motion  is  not  affected  by  a  passing 
shock.  Third,  the  gas  flow  field  is  .decoupled  from  the  par¬ 
ticle  motion.  Fourth,  rocket  motor  motion  does  not  affect 
the  debris  hazard  problem. 

The  first  assumption  is  suggested  by  observations 
of  developing  supersonic  plumes  (see  Appendix  A  on  supersonic 
jets).  Photographs  indicate  that,  at  any  given  time,  a 


developing  plume  can  be  represented  as  a  truncated  version 
of  a  fully  developed  plume  (see,  for  example,  Schmidt,  1974). 

It  follows  that  for  whatever  velocity  a  particle  reaches  in 
a  partially  developed- plume ,  there  will  be  further  accelera¬ 
tion  in  the  fully  developed  plume.  Since  maximum  streamline 
turning  angles  occur  at  the  nozzle, lip,  the  angular  disper¬ 
sion  of  particles  should  be  about  the  same  at  any  time  during 
plume  development. 

The  second  assumption  is  suggested  by  the  fact  that 
a  particle’s  interaction  time  with  shock  waves  emanating  from 
a  small  rocket  is  very  short.  Observations  of  the  motions  of 
ping  pong  balls  in  shock  tubes  -(de  Krasinsky,  1975)  support 
this  assumption. 

The  third  assumption  will  be  valid  when  the  debris 
occupies  an  insignificant  fraction  of  the  plume  volume. 

This  condition  may  be  violated  in  the  vicinity  of  the  nozzle, 
where  debris  loading  of  the  jet  can  be  high.  The  effect  is 
one  of  reducing  jet  momentum  and  deflecting  gas  stream  lines. 

In  any  case,  the  effect  is  expected  to  be  a  small  one  since 
any  alteration  of  flow  field  occurs  only  in  the  nozzle  vicinity. 

The  fourth /assumption  will  be  valid  so  long  as  gas 
speeds  in  the  jet  are  much  greater  than  the  speed  of  the 
motor.  This  condition  will  probably  be  met  in  nearly  all 
cases  —  typically,  a  small  rocket  starts  from  rest'  and  has 
a  terminal  speed  less  than  twenty  percent  of  the  gas  speed. 

With  these  assumptions,  one  can  sefc.^an  upper  limit,  on 
particle  velocity  by  following  the  motions  of  single 
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particles  in  steady  supersonic  jets.  It  is  convenient  to 
transform  equation  (2.1)  from  a  Lagrangian  to  an  Euler ian 
representation  via  the  chain  rule: 

mv^Vvp  =  %CdA  p|vg  -  vp|(vg  -  vp)  (2.2) 

Particle  velocity  v  has  been  transformed  to  spatial  coordin- 
ates  through  the  relation 

t  =  t(x)  (2.3) 

which  is  the  time  t  that  the  particle  occupies  position  x. 

Given  a  jet  flow  field,  equation  (2.1)  can  be  used 
in  conjunction  with  a  skin  penetration  model  (Lewis.,  1978) 
to  predict  a  region  of  debris  hazard.  For  each  class  of 
debris,  one  integrates  equation  (2.2)  to  the  point  where 
(Vpj  drops  below  the  minimum  penetration  velocity  predicted 
by  the  penetration  model.  The  locus  of  all  such  points  de¬ 
fines  the  boundary  of  the  debris  hazard  area. 

Although  the  above  procedure  might  be  useful  as 
a  check  on  more  approximate  procedures,  it  is  not  a  good 
choice  for  a  debris  hazard  model.  As  will  be  shown,  the 
minimum  safe  standoff  distance  is  very  sensitive  to  uncer¬ 
tainties  in  particle  drag  coefficients  and  other  parameters. 

The  computational  effort  required  in  finding  the  flow  field 
and  then  integrating  equation  (2.2)  is  not  justified  by  the 
uncertainty  of  the  results.  Results  of  simpler  models  are 
apt  to  be  just  as  valid. 


I 


The  model  developed  in  the  following  paragraphs 
follows  a  particle  through  three  phases  of  its  motion.  The 
first  phase  concerns  particle  acceleration  in  the  rocket  jet. 
Expressions  will  be  developed  for  estimating  the  particle’s 
maximum  velocity  and  maximum  angular  deviation.  The  second 
phase  concerns  the  particle's  velocity  decay  in  still  air. 

The  third  phase  concerns  the  particle's  impact  with  a  target. 
These  phases  are  shown  schematically  in  figure  (2.1). 

Consider  a  particle  of  mass  m  acted  on  by  a  force 
f(x)  over  a  path  between  two  points  x^  and  X£-  The  speed 
v(x)  of  the  particle  is  changed  only  by  the  component  of  F 
which  is  directed  along  the  path.  That  component  will  be 
designated  F  .  The  particle's  change  of  speed  Av  between 
points  x-^  and  is  then  given  by 


F  2  x 
=  -  / 
m  x* 

*-  o 


1  Fg  (x)  •  dsj  2 


(2.4) 


where  -integration  is  pathwise  and  ds  is  an  element  of 
the  particle's  path  S(x)  between  xQ  and  x-^.  (This  is  equi¬ 
valent  to  integrating  the  left  side  of  equation  (2,2)).  An 
upper  limit  can  be  written  for  Av: 


Av  -  [in1  ^s^max] 


(2.5) 


where  L  is  the  length  of  S  and  |F  I  is  the  maximum  mag- 

°  1  s 'max  ° 

nitude  of  F_. 


For  a  given  path,  Fg  is  equal  to  the  pathwise 
component  of  the  right  hand  side  of  equation  (2.2). 


F„  =  %C  .A.p  i  vrr  -  |  (v„  -  v)  r 

s  d  d  1  g  P  S  p  s 


g 


(2.6) 


where  (v  -  Vp)g  is  the  pathwise  component  of  (v^  -  Vp) . 
Equation  (2.5)  will  continue  to  be  satisfied  if  an  upper 
bound  on  |F  I  is  substituted  for  |F  |  .  From  vector 

S  S  Hi clX 

analysis.,  we  have 


I V  I  I  ( V  -V  )  I  <  ( V  - V  )  *  ( V  “V  ) 

1  g  P1 1 v  g  p's‘-v  g  p'  K  g  p' 


(2.7) 


Substitution  of  (2.6)  and  (2.7)  into  (2.5)  leads  to 


Av 


<  f  ^  C  ,A  p  (v  -v_)  *  (v  -v  )l 

-LM  d  g  P'  g  P  Jmax 


(2.8) 


This  relationship  states  that  the  speed  change  of  a 
particle  is  less  than  it  would  have  been  had  the  particle 
been  acted  on  by  a  constant  force  equal  to  the  global  maximum. 

Relationship  (2.8)  is  much  simplified  when  Vp  is 
assumed  to  be  negligible.  There  is  good  justification  to 
do  this.  In  a  small  rocket,  particles  small  and  light  enough 
to  follow  the  gas  flow  closely  will  have  short  ranges  once 
ejected  from  the  plume.  The  larger  particles,  owing  to  their 
larger  masses,  will  respond  more  sluggishly  to  the  aerodynamic 


forces  of  the  jet  so  that  |v 


|v  |  becomes  the  expected 
S 


<< 


condition.  In  any  event,  the  effect  of  assuming  =  0  is 

i 

generally  to  overestimate  A v  in  equation  (2.4) ,  so  that 
relationship  (2.8)  remains  valid.  With  these  considerations, 
relationship  (2.8)  can  be  written 


4 


CdA  pv? 


(2.9) 


where 


v  •  v 
g  g 


(2.10) 


As  previously  noted,  and  A  can  be  complicated 
functions  of  particle  dynamics  and  aerodynamics.  We  will 
here  suppose  that  constant  values  can  be  chosen  for  each 
such  that  the  inequality  in  (2.9)  is  not  violated.  How 
those  values  might  be  chosen  is  deferred  to  later  discussion. 
Relationship  (2.9)  becomes 


L  CdA 


(2.11) 


The  quantity  pv2  in  relation  (2.11)  is  the  local 

g 

momentum  flux  density  of  the  plume.  If  the  gas  in  the  plume 

2 

undergoes  an  isentropic  expansion,  then  pv  has  a  maximum 

g 

determined  by  the  gas  dynamics.  (The  assumption  of  an  iseri- 
tropic  expansion  may  not  be  valid  if  a  significant  amount 
of  heat  is  added  during  the  expansion.  This  may  occur  in  a 
small  rocket  if  significant  amounts  of  fuel  are  burned  out¬ 
side  the  combustion  chamber.) 


In  an  isentropic  expansion  (Van  Wylen,  p.35$), 
2kRT  /  m  \ 


Vg  =  k 


P  *  P„ 


»  1 

f 

’  lTo  / 


*-(*) 


(2.12) 


(2.13) 


(2.1k) 


where  the  subscript  o  denotes  stagnation  conditions  and 
where  k  is  ratio  of  specific  heats,  R  is  gas  constant,  T 
is  temperature,  and  P  is  pressure.  Substituting  (2.1^1) 
into  (2.12)  and  (2.13)  and  noting  that 

RTp  =  , 

o  o  o 


(2.15) 


one  is  lead  to 


&  te  r ■  t) 


(2.16) 


Equation  (2.16  is  plotted  in  figure  (2.2)  for  k  equal  to 


l.J|,  corresponding  to  air,  and  for  k  equal  to  1.16,  corres¬ 
ponding  to  hot  combustion  gases.  The  figure  shows  variation 

2 

in  the  dimensionless  drag  force, ,  as  the  flow  expands 
from  stagnation  =  1^  to  vacuum  =  0^.  Note  that  the 

peak  values  of  these  curves  differ  little  in  magnitude 


and  occur  at  approximately  the  same  values  of  pressure  ratio. 
This  suggests  tha .•  a  particle  drag  model  depending  on  maximum 


momentum  flux  density  should  be  insensitive  to  uncertainties 


in  k. 


The  rapid  fall  off  in  dimensionless  drag  force 

p 

as  p—  approaches  zero  is  an  important  .consideration  in 
*0 

laser  doppler  velocimetery .  That  technique  assumes  that 

drag  forces  are  large  enough  that  entrained  particles  follow 

the  flow  streamlines  everyv/here .  In  an  underexpanded  super- 

P 

sonic  jet,  where  p—  can  easily  be  as  small  as  0.001, 
particles  may  not  follow  the  flow. 

Differentiating  equation  (2.16)  with  respect  to 

p 

«—  and  setting  the  result  equal  to  zero  leads  to 
*o 


-  00 

o 


(2.17) 


for  the  pressure  ratio  corresponding  the  maximum  momentum 
flux  density.  Substitution  of  (2.17)  into  (2.16)  leads  to 


pv2  =  2P  k 


(2.18) 


for  the  maximum  momentum  flux  density.  Equation  (2.18) 
states  that  the  maximum  momentum  flux,  and  hence  the  maxi¬ 
mum  drag  force  on  a  particle,  is  directly  proportional  to 
the  stagnation  pressure  and  weakly  dependent  on  the  ratio  of 
specific  heats  over  its  normal  range.  One  can  use  equation 
(2.14)  and  the  expression  for  Mach  number  M, 


to  show  that  the  maximum  in  momentum  flux  density  occurs  in 
an  isentropic  expansion  at  a  Mach  number  equal  to  ijT.  This 
is  very  early  in  the  expansion  and  will  probably  occur  within 


the  nozzle. 


Substituting  (2.18)  into  (2.11)  leads  to 


2LC  ,A  P^  I-E 

—A - 2.  k 

m 


(2.20) 


for  an  upper  bound  on  the  velocity  change  of  a  particle 
entrained  in  an  isentropic  expansion  over  a  length  L* 

The  value  of  L  will  depend  on  the  jet*s  geometry. 

We  propose,  tentatively,  that  it  be  set  equal  to  the  length 
A  of  the  bottle  shock,  noting  that,  downstream  of  the  Mach  ' 
disc,  the  momentum  flux  density  is  rapidly  dispersed  by  tur¬ 
bulent  mixing  (Che-Haing,  1969).  This  choice  is  further 
motivated  by  the  fact  that  the  jet  will  scale  on  X  (Che-Haing, 
1969),  so  that  setting  L=A  is  in  error  by  at  worst  a  multi¬ 
plicative  constant.  It  will  be  shown  later  that  safe  standoff 
distance  is  not  critically  dependent  on  L. 

Lewis  (1966)  gives  an  empi.rical  expression  for 
the  wavelength  X  as  a  function  of  exit  Mach  number  Mg  and'  of 
the  ratio  of  exit  pressure  Pg  to  ambient  pressure  P^: 
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U~  =  0 . 69MP 
dE  E 


(2.21) 


where  dg  is  the  nozzle  exit  diameter.  This  equation  has 

P 

been  verified  over  a  wide  range  of  values  for  MP  and  E 

*T ' 

It  might  be  noted  that  Love  (1959)  gives  an  expression 
which  differs  considerably  from  equation  (2.21).  Lewis 
comments  that,  between  investigators,  there  is  frequently 
much  scatter  in  these  data. 

The  exit  Mach  number  can  be  fouijd  from 


“is  (I^T)(1  +  k2L  **) 


(k+l)/2(k-l) 


(2.22) 


where  Ag  is  the  nozzle  exit  area  and  A*  is  the  throat  area. 
The  exit  pressure  ratio  can  be  expressed  as  (Van  Wylen,  1963) 


!e_!o  r 
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(k-l) 


k/(l-k) 


(2,23) 


Once  clear  of  the  bottle  shock  region  (see 
Appendix  A  on  jet  structure),  a  particle  interacts  with  a 
decaying  jet  plume  where  gas  momentum  flux  density  drops 
rapidly  with  increasing  distance  from  the  source.  Those 
particles  of  most  interest  to  the  debris  hazard  problem  are 
unlikely  to  be  further  accelerated  in  this  region  of  rapidly 
dwindling  jet  influence.  Within  a  distance  of  a  few  X  , 
the  more  dangerous  particles  will  be  traveling  much  faster 


than  the  surrounding  gases.  There,  setting  v  =  0  in  equation  (2.2) 

S 

becomes  a  valid  approximation.  (Note  that  if  the  particles  are 
accelerated  to  supersonic  speeds,  they  can  outrun  the  shock  wave 
generated  by  the  starting  jet.  Setting  v„=  0  in  that  case  is  not 
an  approximation,  but  a  statement  of  the  observed  physics .  Schmidt 
(1574)  has  taken  remarkable  photographs  of  such  particles  occurring 
in  the  muzzle  blast  of  an  M-16  rifle.) 

From  equation  (2.2),  the  motion  of  a  particle  through 
still  air  with  no  body  forces  is 


m  v  ”^CdApAv 


(2.24) 


where  p^  is  the  air  density  aid  where  x  has  been  takei  to  be  in 
the  direction  of  particle  notion.  For  m,  Cj,  A  and  p^  constant, 
equation  (2.24)  integrates  to 


CdApA 


(2.25) 


where  vQ  is  the  initial  velocity.  Equation  (2.25)  expresses  the 

distance  required  for  a  particle  to  drop  in  speed  from  vq  to  v- 

Lewis  (1978)  has  shorn  that  a. particle's  probability 

2 

of  penetration  is  a  reonotonicaliy  increasing  function  of 

q.p 

where  m  is  the  particle’s  mass,  v  is  its  velocity,  A  is  its 

area,  and  Q  is  a  parameter  describing  the  target's  material 

properties  (see  Appendix  B) .  It  follows  that  if  the  probability 

.  2 

of  penetration  is  not  to  exceed  some  given  level,  then  :M— 
imst  not  exceed  some  rusiber  n: 
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or 


(2.26) 


(2.27) 


where 


(2.28) 


will  be  termed  the  penetration  parameter. 

Equations  (2.20),  (2.25),  and  (2.27)  can  now  be  com¬ 
bined  to  form  the  basic  model  for  calculating  safe  standoff 
distance  from  a  snail  rocket.  A  safety  criterion  can  be  set  by 
demanding  that  a  particle,  accelerated  to  the  maximum  speed  given 
by  equation  (2.20),  mist  be  slowed  via  drag,  described  in  equation 
(2.25),  below  the  minimum  penetration  speed  given  in  equation  (2.27). 
This  leads  to 


S  > 


m 

2CdApA 


MdPo 

% 


+ 


+ 


L 


(2.29) 


or,  in  dimensionless  expression 


(2.30) 


where  3  is  the  standoff  distance. 

The  angular  dispersion  of  the  debris  could  be  calculated 
by  integrating  equation  (2.2)  over  a  large  number  of  cases.  As 
with  standoff  distance,  however,  uncertainties  in  the  problem  lead 
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us  to  search  for  a  simplified  means  of  establishing  an  upper 
bounds  on  dispersion  angle.  Gas  dynamics  is  helpful  here. 

The  largest  angle  through  which  a  gas  molecule  can 
tum  in  a  supersonic  expansion  is  given  by  the  Prandtl-'fcyer 
turning  angle  (Shapiro,  1953).  Thus,  in  the  plume  .of  arfmder- 
expanded  rocket ,  no  streamline  will  have  an  inclination  to  the 
axis  largo:  than  the  sum  of  the  nozzle  divergence  and  the  Prandtl- 
Keyer  turning  angle  for  conditions  at  the  nozzle  lip.  See  figure 
(2.3) .  In.  the  absence  of  lift  and  body  forces,  an  entrained  par¬ 
ticle  moves  away  from  the  axis  at  no  angle  larger  than  that  of 
the  most  inclined  streamline.  We  vlll  take  this  as  the  upper 
bounds  on  particle  dispersion.  Conditions  under  which  that  "upper 
bounds"  might  be  exceeded  will  be  discussed  in  a  later  section. 

The  turning  angle  9  for  a  Prandtl-loeyer  expansion  is 
given  by  Shapiro  (1953) : 


0  =-atan 


(2.31) 


This  represents  the  angle  through  which  a  jflow  turns  in  expanding 
from  i’fech  number  unity  to  Vsadi  timber  M. 

To  calculate  the  maxintm  angular  streamline  deviation 
6irax  f°r  m  ^ derexpandsd  jet,  the  following  procedure  can  be  fol¬ 
lowed  .  First,  determine  a  k  appropriate  for  the  propellent  gases. 
Next  determine  an  ambient  turning  angle  0^  by  substituting  into 
equation  (2.31),  where 


/  A-//P  \k_1 

„  2  pi  !  cPtt 
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is  the  Mach  ranker  of  a  flow  expaided  fran  clunker  pressure  P 

r  o 

to  arkient  pressure  P* .  Then  determine  a  nozzle  exit  turning 
angle  8g  by  substituting  MU  into  equation  (2.31),  where  MU 
is  the  Hach  ranker  at  the  nozzle  exit  plane,  satisfying 


k-1 

2 


(kfl)/2(k-l) 


(2.33) 


with  A,  being  the  ©it  area  and  A*  the  threat  area  of  the  nozzle. 
£> 

If  8^,  is  the  nozzle  divergence  angle,  then  the  streamline 

deviation  is 


a 

“max 


-  8r 


(2.34) 


Tubs  launched  rockets  require  seme  special  considera¬ 
tions.  A  system  in  which  the  rocket  nozzle's  edit  plane  is  placed 
■  at  cr  near  the  launch  tube's  breech  will  have  essentially  the  same 

debris  hazard  as  the  rocket  alone.  Ehen  the  rocket  is  placed 

J 

higher  up  in  the  launch  tube,  the  launch  tube  can  be  expected  to 
act  as  an  extension  of  the  rocket  nozzle,  Nozzle  od.t  plane  confi¬ 
guration  for  the  debris  hazard  model  is  then  the  launch  tube  breech 
conf igura ticn .  Usually,  this  will  not  nuch  change  the  stadoff 
distance  and  will  narrow  the  dispersion  angle.  Ehen  the  rocket  is 
placed  several  diameters  up  the  launch  tube,  the  dispersion  ray  widen. 
|  This  collates  the  basic  debris  hazard  nodal.  Equation 

•  (2.29),  for  ndrtirajn  safe  standoff  distance,  and  equation  (2.34), 

-  for  soxinim  dispersion,  together  define  a  sector  of  a  circle  to  the 

I  rear  of  a  rocket  mere  a  debris  hazard  nay  err  1st.  There  are  two 


semiempirical  parameters  in  the  model:  Q  ,  describing  the 

Jr 

target  material,  and  L,  describing  the  particle  acceleration 
distance.  Model  usage  is  described  in  Section  3* 
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SECTION  3 

MODEL  USAGE 

This  section  illustrates  use  of  the  debris  hazard 
model  developed  in  Section  2.  The  Viper,  a  shoulder  launched 
antitank  rocket,  is  used  to  aid  in  the  illustration.  A  step 
by  step  procedure  is  presented  for  determining  standoff  dis¬ 
tance  and  dispersion  angle  from  the  pertinent  physical  data. 
Interactive  computer  codes  incorporating  that  procedure  are 
given  in  Appendix  D.  Finally,  a  series  of  dimensionless  plots 
are  shown  which  identify  critical  parameters  in  the  model. 

The  debris  hazard  model  of  this  report  requires 
knowledge  of  eleven  parameters  —  five  to  describe  the 
rocket  motor,  three  to  describe  the  debris,  one  to  describe 
the  target,  and  two  to  describe  the  ambient  air.  Viper 
characteristics  are  summarized  in  Appendix  C,  and  target 
toughness  is  discussed  in  Appendix  B.  (Additionally,  there 
is  a  s.emiempirical  constant  which  relates  debris  acceleration 
length  L  to  jet  primary  wavelength  X.  In  the  present  implemen¬ 
tation  of  the  model,  L  is  taken  to  be  equal  to  X.)  These  para¬ 
meters  are  listed  in  Appendix  C. 
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The  next  few  paragraphs  show  the  use  of  the  para¬ 
meters  of  Table  (C.l)  to  determine  standoff  distance  and  dis¬ 
persion  angle  for  the  Viper.  The  procedure  is  outlined  in 
figures(3.1)  and  (3.2). 

The  first  calculation  is  the  determination  of  jet 
primary  wavelength  X  from  equation  (2.21).  This  equation  re¬ 
quires  the  exit  Mach  number  Mg,  which  can  be  determined  from 
equation  (2.22)  (see  Appendix  E),  and  the  exit  plane  pressure 
Pg,  which  can  be  determined  from  equation  (2.23).  (Equations 
(2.22)  and  (2.23)  are  represented  in  figures  (3*3)  and  (3-4)  •) 
For  Viper,  =  2.07,  Pg  =  8.03  x  10 6  N/m2,  and  X  =  0.85  m. 
Substitution  of  these  values  and  of  values  from  Table  (C..1) 
into  equation  (2.29)  gives  a  minimum  safe  standoff  distance 
S  =  33m  for  Viper’s  detente  fingers  striking  people  wearing 
summer  weight  uniforms . 

Calculation  of  dispersion  angle  is  as  follows. 
Equation  (2.31):.  using  Mg  determined  in  the  standoff  calcula¬ 
tion,  gives  6-,  =  35°.  Equations  (2.32)  and  (2.31)  give 
Ma  =  *J.27  and  0^  =  102°.  (Equations  (2.32)  and  (2.31)  are 
represented  in  figures  (3-15)  and  (3*6).)  Substitution  of 
0E,  0A  and  0^  (from  Table  C.l)  into  equation  (2.3^)  gives 

6max  ‘  ,8°; 

In  summary,  the  model  predicts  that  a  hazardous 
region  exists  behind  the  rocket  to  a  distance  of  33  meters 
from  the  nozzle  and  to  78  degrees  off  axis.  In  a  complete 
analysis,  the  standoff  distance  would  be  calculated  for  each 
type  of  debris.  The  procedures  outlined  above  have  been  in¬ 
corporated  into  the  computer  codes  documented  in  Appendix  D. 


In  one  of  Viper's  proposed  configurations,  the 
launch  tube  is  extended  rearward  several  centimeters  beyond 
the  rocket  nozzle  exit  plane.  In  that  case,  it  is  approp¬ 
riate  to  regard  the  launch  tube  as  an  extension  of  the  nozzle 
The  nozzle  exit  diameter  d£  should  then  be  set  to  the  launch 
tube  inside  diameter  of  0.0,7-93  m,  and  the  nozzle  divergence 
angle  should  be  set  to  0°.  Calculating  debris  hazard  as 
before,  one  finds  for  andoff  distance,  S  «  3*i  m,  and  for 
off  axis  divergence  angle,  6max  =  52°.  Thus*  a  short  launch 
tube  extension  has  little  effect  on  standoff  distance,  but 

reduces  the  dispersion  angle. 

It  might  be  noted  that,  in  the  above  example,  S 
is  much  greater  than  L.  This  will  generally  be  the  case 
if  there  are  particles  accelerated  to  speeds  well  in  excess 
of  that  necessary  to  cause  damage.  The  error  introduced  by 
dropping  the  last  term  in  equations  (2.29)  and  (2.30)  is, 
therefore,  apt  to  be  small.  With  these  considerations,  one 
can  neglect  that  term  and  express  minimum  standoff  distance 
by  rewriting  equation  (2.30): 
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Equation  (3.1)  is  plotted  in  figures  (3.7)  and  (3-8) .for  , 
several  drag  coefficients  and  for  specific  heat  ratios  of 
1,1b  and  1.67.  The  variation  of  equation  (3-1)  with  lc  over 
the  range  l,15<k<1.6?  is  slight— see  figure  (3-7).  Note, 
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however,  that  there  is  an  implicit  variation  with  k 
contained  in  L. 

Figure  (3*7)  covers  a  parameter  range  likely 

for  small  rockets  and  unarmored  personnel.  Two  points 

are  noteworthy  here.  First,  the  standoff  distance  is 

LP 

only  weakly  a  function  of  L  for  —x—  greater  than  about 

QP 

20.  This  means  that  the  estimated  standoff  distance  is 
not  much  affected  by  errors  in  estimation  of  L.  Such 
behavior  is  generally  desirable  for  semiempically  deter¬ 
mined  parameters.  Second,  the  standoff  distance  is  very 
strongly  a  function  of  drag  coefficient.  Since  drag  co¬ 
efficients  can  vary  sharply  with  particle  geometry  and 
dynamics  (Sadeh,  1975),  this  causes  a  serious  and  unavoid¬ 
able  uncertainty  in  the  debris  hazard  model.  Specifica¬ 
tions  of  drag  coefficients  are  likely  to  be  the  overriding 
source  of  undertainty  in  any  debris  hazard  model.  There¬ 
fore,  further  sophistication  in  the  jet  fluid  dynamic 
model. is  likely  to  be  unproductive.  The  reader  is  refer¬ 
red'  to  the  discussion  in  Section  4. 


Figure  (3*8),  for  small  values  of 
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,  covers  a 
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range  that  might  be  pertinent  to  armored  personnel.  The 

LP. 

standoff  distance  for  this  range  varies  rapidly  with 


o 


and  chaotically  with  drag  coefficient.  Clearly,  a  large 
margin  of  safety  should  be  applied  for  positions  in  this 
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range. 
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Figure  (3*9)  illustrates  the  quantity  0  =  (0.-6=) 

,  0  ii 

as  a  function  of  ambient  pressure  ratio  P^/P^  ,  exit  to 

* 

throat  nozzle  diameter  ratio  d^/d  ^  and  three  values  of 
specific  heat  ratio.  The  angle  0Q  represents  the  maximum 
off  axis  divergence  of  fluid  stream  iines  for  a  straight 
supersonic  nozzle.  These  figures  can  be  used  in  lieu  of 
evaluating  equation  (2.31).  Evidently,  0Q  is  a  strong 

g 

function  of  all  three  of  its  arguments,  Pq/Pa,  dg/d  * 


and  k. 


Helium,  with  a  mass  density  near  that  of  typical 


combustion  products,  is  occasionally  used  in  modeling 
chemical  rocket  jets.  But,  helium  has  a  specific  heat 
ratio  of  1.67-  The  strong  dependence  of  6q  on  k  argues 
that  caution  be  used  in  deriving  quantitative  conclusions 
from  Such  experiments. 
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Figure  (3.1).  Outline  for  Calculation  of  Standoff. Distance. 


For  standoff  distance  S,  use  equation  (2.29): 
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2LC,P0 

ln  ~Q^  +i^Llnk  +L 


A)  For  m,  Cd>  A,  p^,  PQ,  and  ka  see  Table  .(3*1) 


B)  For  skin  penetration  parameter  Q  ,  see  Appendix  B 

P 


C)  For  acceleration  length  L  use  equation  (2.21) 


/kP  \h 

L  =  X  =  0.69  Hgdg 


1)  For  dj,,  k  and  P.,  see  Table  (3-1) 

£»  n 


2)  For  exit  Mach  number  M£,  use  equation  (2.22) 


(^S.  V  =  if  /JlNA  +  £=1 1/2  ) 
\ ds  /  me  lk+i A  2  r*E  / 


(k+l) 
2Tk =1) 


a)  For  d*,  see  Table  (3.1) 

b)  See  Appendix  E  for  inversion  of  equation (2. 22) 

3)  For  exit  plane  pressure  P-,  use  equation  (2.23) 


Pp  =  P„ 
E  o 
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Figure  (3*2).  Outline  for  Calculation  of  Dispersion  Angle, 


For  dispersion  half  angle  0  use  equation  (2.3*0 

0  __  =  0fl  -  0„  +  0., 

max  A  &  N 


A)  For  0^,  see  Table  (3.1) 

B)  For.  0E,  substitute  Mg  into  equation  (2.31) 

0E  -  -atan(H|  -  ljH  +  0f  atan  ^(m|  -  1  j) 

1)  For  k,  see  Table  (3*1) 

2)  For  exit  Mach  number  Mg,  use  equation  (2.22) 


k-ll,2\jk+1/2(k"1) 
2  l  iE 


(^=1  /_J_V3 

\d«;  Mg  yk+1 


a)  For  dg  and  d*  ,  see  Table  (3*1) 

b)  See  Appendix  E  for  inversion  of 
equation  (2.22) 

C)  For  0^.,  substitute  into  equation  (2.31) 

8J i*  -atan(MA  -  l)  S  -(gf)’"  ^ 

1)  For  use  equation  (2.32) 


a)  rur  and  P. ,  see  Table  (3*D 
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adiabatic  expansion 


abafcic  expansion. 


turning  angle,  versus  Mach  number.,, 
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3 f  constant  nozzle  lip  turning  angles  plotted 


MODEL  RESTRICTIONS 


Tills  section  discusses  conditions  in  which  the 
model  of  this  paper  may  underestimate  the  extent  of  "'‘he 
debris- hazard  area.  Such  conditions  can  arise  when  the 
input  parameters  of  the  model  do  not  adequately  describe 
the  particles  or  the  jet.  The  following  paragraphs  will 
discuss  the  effects  of  particle  shape,  obstructions,  and 
external  burning. 

It  is  assumed  in  the  model  that  a  particle  of 
debris  is  .completely  specified  by  its  mass,  projected  area, 
and  drag  coefficient.  This  description  may  be  inadequate’ 
for  some  classes  of  particle  shapes.  Particles  with  one 
very  long  or  very  short  dimension  (such  as  a  length  of  wire 
or  a  thin  disc)  can  develop  lift  and  drag  forces  far  more 
complex  than  those  of  the  simple  model  assumed  here.  Par¬ 
ticles  with  protruding  points  or  edges  may  have  penetrating 
powers  well  in  excess  of  that  suggested  by  Lewis*  (1978) 


model. 

As  noted  in  the  proceeding  section,  drag  co¬ 
efficient  is  a  critical  parameter  in  the  model.  Unfortun¬ 
ately,  It  is  a  parameter  difficult  to  estimate  for  particle 
of  Irregular  shape.  A  problem  can  arise  if  a  particle  pre¬ 
sents  a  large  area,  high  drag  coefficient  facet  to  the  jet 
flow,  but  reorients  itself  and  presents  a  small  area,  small 


drag  coefficient  facet  for  its  motion  through  the  still 
air.  Such  a  particle  can  have  a  greater  range  than  a 
similar  one  described  by  a  single  drag  coefficient. 
Equation  (2.24),  giving  the  minimum  safe  standoff  dis¬ 
tance,  can  be  rederived  keeping  drag  coefficients  and 
projected  areas  separate  for  different  phases  of  particle 
flight.  In  that  case. 
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where  C - .  and  A.  are  drag  coefficient  and  projected  are 
for  motion  in  the  jet,  C .  and  A_  for  motion  in  the  still 
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air,  and  A  is  the  projected  area  pertinent  to  particle 
* 

impact.  Equation  (4.1)  is  probably  of  academic  concern 
only,  since,  for  a  given  class  of  particle,  the  semi- 
empirical  constant  L  can  absorb  the  differences  in  drag 
coefficients  and  projected  areas. 

Lift  forces  act  perpendicularly  to  a  particle's 
relative  wind,  and  can  serve  to  increase  the  angular  dis¬ 
persion  of  debris  from  a  rocket  backblast.  (Baker  (1978) 
discusses  a  lift  induced  range  increase  for  debris  from 
explosions.  This  can  occur  if  lift  forces  are  so  oriented 
as  tc  flatten  a  particle's  trajectory.  Small  particles, 
as  might  be  expected  in  a  small  rocket  backblast,  will 
slow  to  safe  soeeds  before  the  range  lengthening  effects 


of  lift  can  be  felt.  In  any  case,  the  model  of  this  paper 
assumes  the  worst  case  condition  of  no  trajectory  curvature.) 
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The  angular  dispersion  model  of  this  paper  is  conservative 
in  the  sense  that  it  will  overestimate  particle  dispersion 
in  the  absence  of  lift.  Careful  consideration  of  fluid 
mechanics  will  reveal  that  a  particle  moving  at  the  maxi¬ 
mum  predicted  angular  deflection  will  have  received  zero 
momentum  in  that  direction.  This  should  insure  an  adequate 
margin  of  safety  for  irregularly  shaped  particles.  A  thin, 
flat  object,  however,  can  literally  fly  out  of  the  pre¬ 
dicted  region  of  hazard  in  a  highly  unpredictable  fashion. 

If  such  objects — nozzle  diaphrams  or  end  caps,  for  example — 
are  expected  to  be  present  in  the  debris,  they  should  be 
considered  to  present  a  360°  hazard  until  testing  proves 
otherwise.  The  reader  is  reminded  that  a  boomerang  can 
readily  negotiate  a  720°  turn  while  maintaining  enough 
momentum  to  inflict  appreciable  damage. 

It  is  assumed  in  the  model  that  the  region  behind 
the  rocket  is  clear.  Obstructions  In  the  predicted  region  cf 
debris  hazard  can  deflect  debris  either  directly  through 
ricochet  or  indirectly  through  jet  deflection  or  spalling. 
Ricochet  can  be  expected  to  increase  particle  dispersion, 
possibly  to  the  point  of  endangering  the  gunner.  Frangible 
objects  within  a  few  jet  wavelengths  of  the  rocket  nozzle 
may  not  only  deflect  debris,  but  may  break  up  ana  themselves 
become  part  of  the  debris  hazard.  Again,  these  effects  are 
situation  dependent  and  likely  to  be  unpredictable.- 


Particle  range  is  unlikely  to  be  much  affected.  Par¬ 
ticle  dispersion  can  be  expanded  to  360°  if  there  are 
solid  obstructions  in  the  near  field. 

It  has  been  assumed  that  the  rocket’s  jet  is 
an  adiabatic  process,  that  is,  that  any  combustion  occurs 
within  the  rocket.  The  jet  geometry  may  be  altered  if 
afterburning,  or  combustion  external  to  the  rocket#  takes 
place.  An  extreme  example  of  this  is  the  after-burn  ex¬ 
plosion#  which  may  occur  when  substantial  amounts  of  un- 
bur-ned  propellant  accumulate  in  the  plume  before  igniting. 
Afterburning  in  tactical  rockets  has  been  studied  for 
many  years  ( JANNAF,  1976).  It  is  a  highly  complex  phe¬ 
nomenon!,  depending  on  the  interaction  of  propellant-  che¬ 
mistry,  ambient  conditions,  and  fluid  dynamics.  Kuch  is 
as  yet  poorly  understood. 

In  general,  one  can  expect  afterburning  to  in¬ 
crease  both  the  range  and  dispersion  of  backblast  debris. 
The  increase  in  range  is  likely  to  be  small  for  two  reasons 
First,  energy  release  in  an  afterburn  is  probably  only  a 
fraction  of  the  energy  release  in  a  rocket’s  combustion 
chamber.  Second,  debris  range  for  small  rockets  is  relat¬ 
ively  insensitive  to  changes  jet  plume  characteristics. 

This  is  true  fox’  that  same  reason  that  standoff  distance 
varies  only  slowly  with  L  over  the  range  of  interest  (see 
figure  (3-7)).  Dispersion,  on  the  other  hand,  may  be 
greatly  increased,  particularly  in  the  event  of  an  after¬ 
burn  explosion.  An  after-burn  explosion  can  increase 
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dispersion  in  at  least  two  ways:  it  can  deflect  particles 
already  accelerated  in  the  main  jet,  and  it  can  itself  ac¬ 
celerate  particles.  The  nature  of  these  explosions  is  not 
understood,  either  analytically  or  empirically,  to  an  extent 
that  would  allow  reliable  estimates  to  be  made  for  dispersion 
from  a  given  rocket.  If  afterburning  is  expected  to  occur 
in  a  given  system,  tests  should  be  designed  to  detect  debris 
with  anomolously  large  dispersion  angles.  In  the  Viper  test 
program,  for  example,  particles  of  des leant  have  been  caught 
on  collection  panels  placed  at  the  gunner’s  position  (Chipser). 
Presumably,  those  particles  were  blown  forward  by  energetic 
events  occurring  behind  the  launch  tube. 

While  the  debris  hazard  model  has  not  yet  been 
verified  quantitatively  against  experiment,  some  qualita¬ 
tive  observations  were  made  during  test  firings  of  LAW  and 
Viper  at  the  Redstone  Arsenal  during  1980.  Both  of  these 
weapons  are  shoulder-fired,  tube  launched  antitank  rockets 
of  the  sort  for  which  the  model  Is  designed.  Viper  Is 


described  in  more  detail  in  Appendix  C. 

iijg  pi) to  iPti^oducG  ^  j_o  P'rsin  s  tGcl 

cubes  into  the  backblast  of  one  LAW  firing-  Those  cubes . 


simiiai 


;o  those  used  by  Lewis  (1978),  were  glued  to  veil- 


paper  stretched  across  the  breech  of  the  launch  tube. 
Witness  panels,  consisting  of  two  Inches  of  foamed  plastic 
backed  by  half  an  inch  of  plywood,  were  placed  at  ten 


meters  and  twenty  meters  directly  behind  she  LAW.  A 
streak  camera  was  set  up  to  record  particle  velocities 
approximately  one  half  meter  behind  the  launch  tube 
breech.  Three  cubes  were  recovered  after  firing — two 
from  the  plywood  of  the  ten  meter  panel  and  one  from  the 
plywood  of  the  twenty  meter  panel.  There  was  also  a 
hole  through  the  ten  meter  panel  'where  an  object,  presum¬ 
ably  a  cube,  had  penetrated.  A  number  of  dense  streaks 
were  recorded  by  the  camera,  corresponding  to  particles 
moving  to  the  rear  between  150  ana  h‘50  meters  per-  second. 
There  is  no  way  to  determine  whether  any  of  the  streaks 
were  caused  by  cubes.  Additionally,  there  was  one  streak 
which  would  have  been  caused  by  a  particle  moving  forward 
from  some  place  to  the  rear  of  the  LAW  with  a  velocity  of 
approximately  200  meters  per  second.  We  have  speculated 
that  this  may  have  been  a  particle  propelled  by  an  after¬ 
burn  explosion.  The  witness  panels  have  never  been  cali¬ 
brated  to  correlate  particle  penetration  to  particle  velo 
city  or  to  skin  penetration.  It  was  noted,  however,  that 
the  steel  cubes  had  more  penetrating  power  than  other, 

om  the  rocket  motor.  This 


n  agreement  with  tne  oredre 


wii*; 


Again,  these  observations  are  qualitative,  since  the 
experiments  were  uncontrolled  and  uncalibrated  with  regard 
to  debris  hazard.  The  debris  consisted  of  irregular  plas¬ 
tic  "detente  fingers"  about  a  centimeter  long,  pieces  of 
wiring  and  ignitors,  fragments  of  the  plastic  nozzle  closure 
ranging;  in  size  to  the  full  disc,  styrofoam  throat  plugs, 
and  so  forth.  There  was  a  great  deal  of  variation  in  the 
rocket  configuration  from  test  to  test . 

Witness  panels  as  described  above  were  placed 
at  ten  and  twenty  meters  on  some  runs  and  at  the  gunner's 
position  on  at  least  one  set  of  runs.  As  might  be  expected, 
particle  penetration  on  the  ten  meter  panel  was  generally 
greater  than  on  the  twenty  metei  panel.  There  have  been 
some  measurements  made  of  penetration  depth  and  particle 
characteristics,  which  may  be  useful  if  someone  ever 
calibrates  the  panels.  While  the  gunner's  position  panel 
was  intended  to  pick  up  debris  thrown  backward  as  the 
rocket  exited  tne  launch  tube,  it  also  picked  up  desicant 
spheres  (about  one  mm  diameter)  on  the  rear  facing  side 
(Chipser,  1980).  Again,  we  speculate  that  this  was  debris 
blown  forward  by  an  afterburn  explosion. 

Debris  dispersion  „an  be  estimated  from  the 
distribution  of  .launch  site  ground  debris.  After  a  number 
of  firings  of  both  LAV/  and  Viper,  there  was  a  good  deal  of 
debris  found  along  lines  extending  backwards  from  the 
launch  tube  breech  at  about  forty-five  degrees  off  axis. 


Very  .little  debris  was 


found  more  than  sixty  degrees  off  axis. 


While  these  observations  are  substantially  in  agreement 
with  the  model  predictions  (see  Section  3),  we  stress 
their  qualitative  nature.  The  firings  were  uncontrolled 
experiments  in  this  regard,  and  personal  activity  in  the 
launch  site  vicinity  would  have  redistributed  the  debris 
to  some  extent.  John  Chipser  ( 19 80 )  of  the  Human  Engineer¬ 
ing  Laboratory  has  noted  that  the  witness  panels  collect 
a  surprisingly  small  amount  of  debris  at  their  standard  on 
axis  locations.  He  has  recommended  that  the  panels  be 
moved  off  axis,  noting  that  debris  concentration  is  heav¬ 
iest  along  the  aforementioned  forty-five  degree  lines. 

Such  experiments  are  as  yet  pending. 
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APPENDIX  A 


JET  PLUME  STRUCTURE 


This  appendix  presents  an  overview  of  jet  plume 
structure  pertinent  to  the  debris  hazard  problem.  The 
descriptions  will  be  generally  qualitative  with  references 
given  to  quantitative  analyses  where  needed.  Of  primary 
importance  for  the  debris  hazard  model  is  the  extent  of 
the  bottle  shock  region.  An  empirical  relation  will  be 
given  for  this.  Additionally,  some  basic  properties  of 
a  tube  confined  supersonic  flow  will  be  discussed. 

Underexpanded  supersonic  jets  have  complex 
structures  that  have  so  far  defied  any  simple  explana¬ 
tions.  The  method  of  characteristics  (Shapiro,  1953) 
and  an  assortment  of  hydrocodes  (JAKNAF ,  1976)  have  all 
been  used  in  the  last  twenty  or  thirty  years  to  map  out 
velocity  and  pressure  flow1'  fields  in  such  jets.  Con¬ 
currently,  jets  have  been  probed  and  photographed  over 
a  wide  range  of  physical  conditions.  The  enormity  of 
the  problem  is  underscored  by  the  fact  that,  under  tre¬ 
mendous  impetous  from  the  missile  development  programs, 
theory  and  experiment  are  only  beg5.nning  to  come  into 
good  agreement. 


Figure  (Al)  illustrates  the  bottle  shock  region 
of  a  steady,  underexpanded,  supersonic  jet.  At  the  exit 
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plane,  the  jet  has  a  Mach  number*  of  i-U  and  a  pres- 

Hi 

sure  of  P„.  It  expands  into  still  air*  with  a  smaller 

ih 

pressure  of  P^. 

For  now,  assume  the  flow  is  of  a  perfect  gas 
with  no  heat  addition.  In  that  case,  there  will  be  a 
Prandtl-Meyer  expansion  fan  attached  to  the  nozzle  lip. 
This  is  illustrated  in  the  upper  portion  of  figure  (Al). 
The  lines  of  the  fan  represent  characteristics  across 
which  fluid  thermodynamic  properties  change  by  some 
fixed  factor.  As  it  crosses  these  characteristics,  the 
gas  turns  away  from  the  jet  axis  and  drops  in  pressure, 
so  that  a  diverging  source  type  flow  is  developed  in  the 
vicinity  of  the  nozzle.  The  innermost  characteristic  can 
be  interpreted  as  the  disturbance  from  the  nozzle  lip 
propagating  into  the  jet.  This  disturbance  propagates 
inward  at  the  Mach  angle . 

It  is  the  diverging  nature  of  the  flow  that 
eventually  leads  to  computational  problems.  Gas  in  a 
diverging  supersonic  flow  must  continuously  drop  in 
pressure,  but  the  outer  boundary  of  the  jet  is  held  at 
atmospheric  pressure.  The  apparent  paradox  is  resolved 
in  the  flow  when  the  expansion  fan  reflects  from  the 
surface  of  specified  pressure  which  marks  the  jet’s 
boundary.  (For  inviscid  fluids,  that  boundary  is  a 
vortex  sheet,  that  is,  a  velocity  discontinuity.)  The 
reflection  changes  the  sense  of  the  expansion  waves  in 
that  tney  become  compression  waves  with  opposite  their 
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original  turning  angles.  Problems  arise  when  this 
system  of  compression  waves  coalesce  to  form  an  oblique 
shock  imbedded  in  the  plume.  This  shock,  because  of 
its  shape,  is  referred  to  as  a  bottle  shock.  Gas 
within  the  bottle  shock  remains  in  a  divergent  compres¬ 
sible  source  flow.  At  the  shock,  it  turns  sharply  down¬ 
stream  and  forms  a  lower  speed,  though  still  supersonic, 
sheath  b_ tween  the  bottle  shock  and  the  ambient  air. 

Further  complications  arise  downstream,  since 
the  barrel  shock  bends  around  and  reflects  from  the  jet 
axis.  For  overpressure  ratios  greater  than  about  2, 
this  is  generally  a  Mach  reflection  (Love,  1959),  with 
a  Mach  stem  extending  from  the  barrel  shock  towards  the 
axis  to  form  a  so-called  Mach  disc  (or  Rieman  wave  in 
plane  jets).  As  the  Mach  disc  can  be  of  enormous  strength 
and  so  drastically  alter  the  fluid  motion,  its  location 
is  of  great  importance  in  describing  the  jet  --  many 
features  of  the  jet  scale  on  the  primary  wavelength, 

3>e*the  distance  to  the  Mach  disc. 

The  primary  wavelength  Can  be  calculated 
(sometimes!)  with  hydrocodes  or  with  any  of  a  number  of 
empirical  formulas.  It  is  a  function  of  Mach  number  and 
ratio  of  specific  heats,  and,  to  a  lesser  extent,  of 
nozzle  divergence,  combustion  chemistry,  and  particulate 
loading.  Lewis  (1966)  gives  the  relationship 


.69 

& 


LpiY 

>’a  / 


(1  -h  0.197  Mr 


1.45  0.6 


(A.  1) 


where  d  is  nozzle  diameter,  Mg  is  Mach  number  at  nozzle 
exit,  Pg  is  pressure  at  nozzle  exit,  is  ambient  pres¬ 
sure,  and  $  is  fractional  particulate  mass  density.  Lewis 
notes  that  this  fits  his  data  to  within  five  percent,  but 
that  there  is  considerable  scatter  between  data  from 
different  investigations. 

Beyond  the  Mach  disc,  the  flow  takes  on  an 
annular  form.  Gas  passing  through  the  Mach  disc  is 
shocked  subsonic,  leaving  a  subsonic  core  flow  surrounded 
by  a  supersonic  sheath.  Virtually  all  the  momentum  flux 
is  carried  in  this  sheath.  Viscous  and  turbulent  pro¬ 
cesses  transport  momentum  from  the  sheath  and  into  the 
core  and  the  ambient  atmosphere.  Che-Haing  (196S)  has 
treated  these  processes  seraiempiricaliy ,  but  the  results 
are  difficult  to  use.  As  a  general  rule,  velocities  in 
the  jet  beyond  the  Mach  disc  decay  roughly  as  the  in¬ 
verse  of  the  downstream  position  (JAMN.4P,  1976).  This 
occurs  as  momentum  is  transferred,  primarily  through  tur¬ 
bulent  mixing,  to  the  ambient  air. 

An  impulsive  starting  flow  for-  a  supersonic 
jet  is  a  highly  complex  phenomenon.  In  general,  there 
will  be  a  supersonic  region  extending  from  the  nozzle. 
Schlieren  photographs  (Schmidt,  1974)  show  this  region 
to  resemble  closely  a  truncation  of  the  corresponding 
fully  developed  flow.  At  the  downstream  end  of  this 
region  there  will  be  an  expanding  recirculation  zone 
which  entrains  ambient  air  and  interfaces  jet  momentum 


flux  to  the  atmosphere 


There  is  frequently  a  well 


defined  vortex  ring  generated  by  this  process.  Finally 
the  expanding  gases  will  generate  compression  waves  in 
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the  ambient  air  which  can  be  expected  to  form  one  or 


more  spherical  shocks. 

A  tube  confined  supersonic  jet  can  be  profoundl 
different  from.'  a  free  jet.  Scaling  in  such  a  flow  de¬ 
pends  not  only  on  rocket  exit  conditions,  but  also  on 
the  relative  size  of  the  tube,  on  the  distance  to  the 
tube  exit,  and  on  the  characteristics  of  the  tube  wall. 


In  tube  launched  weapons,  the  rocket  nozzle  exit  diameter 
is  probably  not  much  smaller  than  the  inside  diameter  of 
the  tube.  Three  cases  can  be  delineated  according  to  the 
distance  the  rocket  is  placed  from  the  end  of  the  tube. 

For  a  nozzle  position  very  close  to  the  tube 
exit,  the  jet  is  virtually  unaffected.  If  a  line  drawn 
downstream  from  the  nozzle  lip  at  the  maximum  stream¬ 
line  divergence  angle  (see  Section  2)  doer  not  intercept 
the  tube,  then  this  condition  is  satisfied. 

For  nozzle  positions  up  to  a  diameter  or  so 


fruiii  the  tube  exit,  the  tube  can  act  much  as  an  extension 
of  the  nozs3e.  ’’Exit”  conditions  can  be  referred  to  the 
tube  exit  as  opposed  to  the  rocket  nozzle  exit.  Debris 
dispersion  angle  is  likely  to  be  reduced  since  the  over¬ 
pressure  ratio  measured  at  the  tube  exit  is  smaller  than 
the  overpressure  ratio  measured  at  the  nozzle  exit  of  the 


rocnet . 


unconfined 


Poi*  nozzle  positions  further  than  a  diameter 
from  the  tube  end,  there  can  be  a  substantial  boundary  layer 
build  up  and  internal  shocking  inside  the  tube.  These 
phenomena  depend  critically  on  tube  roughness  and  on  the 
nature  of  the  jet  itself.  Conditions  at  the  tube  exit  are 
best  found  empirically  in  this  case.  Once  those  conditions 
are  known,  the  nature  of  the  external  plume  can  be  estimated 
by  previously  outlined  procedures.  It  is  important  to  the 
debris  hazard  problem  that  the  pressure  at  the  tube  exit 
in  this  case  can  be  larger  than  the  pressure  at  the  nozzle 
exit  of  the  unconfined  rocket.  Contrary  to  intuition,  extend¬ 
ing  the  breech  of  a  rocket  launch  tube  has  the  potential  of 
increasing  debris  dispersion  angle. 
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APPENDIX  B 


PENETRATION  OF  SKIN  BY  BALLISTIC  PARTICLES 
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Lewis  (1978)  has  developed  a  semiempirical  model 
for  estimating  the  probability  of  skin  penetration  by  high 
speed  debris.  A  series  of  experiments,  in  which  objects 
were  thrown  against  simulated  -skin/clothing  combinations, 
was  used  to  calibrate  the  model.  The  objects  included 
tungsten  and  steel  cubes,  wooden  cylinders,  and  gravel; 
their  masses  were  typically  between  a  quarter  gram  and 
four  grams .  Three  skin/clothing  combinations  were  investi¬ 
gated,  simulating  bare  skin,  summer  weight  uniform,  and 
winter  weight  uniform.  Lewis 1  report  contains  extensive 
tables  of  the  experiment  results.  The  model  is  presented 
here  in  a  somewhat  modified  form. 

It  can  be  reasoned  that  skin  penetration  must  be 
a  function  of  particle  dynamics  and  geometry  and  of  target 
material  properties.  For  geometrically  similar  particles 
and  homogeneous  target  material,  a  model  must  certainly 
include  a  particle  characteristic  length  D,  particle  mass  m 
particle  velocity  v,  and  an  areal  toughness  T.  (Areal 
toughness  is  defined  as  the  amount  of  mechanical  energy 
that  can  be  absorbed  per  unit  area  of  a  thin  material  befor 
fracture  occurs.  It  has  units  of  energy  per  length  squared 
These  parameters  can  be  combined  to  form  a  dimensionless 


n  is  evidently  a  measure  of  fch’e  kinetic  energy 
delivered  per  unit  area  of  target  divided  by  the  target 
areal  toughness-  We  expect  that  for  II  large  enough,  the 
mechanical  energy  deposited  in  the  target  exceeds  the 
target's  capacity  to  absorb  energy  without  fracture. 

In  general,  the  particle  is  irregular  and  tumbling, 
and  the  target  is  inhomogeneous,  anisotropic,  and  rate 
dependent.  Neglecting  these  factors  will  result  in  an- 
apparent  scatter  of  the  n*enetration  data.  For  example,  if 
the  particle  is  a  tumbling  cube,  there  may  well  be  a  range 


for  IT  in  which  penetration  occurs  when  the  cube  impacts 
with  a  corner,  but  does  not  occur  ’when  the  cube  impacts 


with  a  flat.  A  probabilistic  model  is  then  appropriate: 


p  {skin  penetration}  =  f(JI) 


(32) 


*{E}  is  the  probability  density  function  for  occur- 


wiicrc  u 


a  nee  or  event;  n  - 


Lewis  has  fitted  S  curves  fee  his  data  using 
the  Walker-Duncan  method  (Walker,  1967).  He  uses  the  par- 
tide's  projected  area  A  instead  of  D  .  These  curves 
are  replotted  in  figure  (Bl).  The  ordinate  is  the  natural 

.  a 

£ 

logarithm  of  the  numerical  value  of  in  SI  units.  (Note 
that  Lewis  uses  a  system  of  mixed  units.)  As  a  result,  a 
different  curve  is  obtained  for  each  value  of  T,  that,  is, 
for  each  skin/clothing  configuration.  The  50%  probability 
level  Is  summarised  in  table  (Bl). 


numerical 


APPENDIX  C 


VIPER  CHARACTERISTICS 

The  Viper  is  a  shoulder  launched  light  anti¬ 
tank  rocket  under  development  by  General  Dynamics.  The 
motor,  using  a  high  performance  boron  based  propellent, 
burns  out  within  the  launch  tube.  Data  given  below 
apply  to  a  configuration  tested  in  the  spring  of  1980, 
and  may  not  be  representative  of  current  design.  They 
are  pr>_  ?nted  to  illustrate  the  sort  of  data  necessary 
to  use  the  debris  hazard  model. 


TABLE  C . 1 


VIPER  CHARACTERISTICS 


PARAMETER 


Rocket  Motor: 


SYMBOL 


VALUE  FOR  VIPER 


Nozzle  exit  diameter,  .dg . 0 . 0 2 iJin 

Nozzle  throat  diameter. d* . 0.06l6m 

Nozzle  divergence  .  .  .GM . 11° 

angle 

Chamber  pressure.  .  .  .P . 6.8x10'  N/m£ 

c 


Specific  heat  ratio  .  .k . 1.16 

for  combustion 
products 


Debris  (plastic  detente  fingers  used  to  secure  rocket 

in  launch  tube) : 

2 

Projected  area . A? . 10  m 

_!( 

Mass . m  .  .  . . 7x10  k£ 

.  Drag  coefficient.  .  .  .  C^ . 0.5 


1’arget  (summer  weight  uniform  -  two  layers  of  clothing): 

5  2 

Skin  penetration.  .  .  .Q . .  .  .  3-^xlO  kg/sec 

parameter  ^ 


Ambient  air  (sea  level,  25"C): 


Mass  density . p^.  .  .  . 

Pressure . P . 

A 


.  1.29  kg/m' 

E,  O 

.  10"  n/m 


Bin  . nniiiit . . . 
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COMPUTER  AND  CALCULATOR  CODES 


Two  codes  are  given  here  which  implement 
the  debris  hazard  model  as  developed  in  Section  2. 
The  first  code  is  written  in  floating  point  BASIC 
and  should  run  on  any  machine  supporting  that  langu¬ 
age.  The  second  code  is  written  for  the  Hewlett- 
Packard  ^lC  programmable  calculator  and  is  specific 
to  that  device.  The  algorithms  used  in  these  codes 
are  discussed  in  Program  Notes  sections,  and  can  be 
adapted  to  other  languages. 


DAO KB LAST  DEBRIS  HAZARD  CODE 


£■  n 

i// 


IN  BASIC 


PROGRAM  NOTES: 


This  code  Is  documented  by  remarks  con¬ 
tained  -in  the  listing.  Note  that  the  listing  appear¬ 
ing  in  this  appendix  io  for  SI  units.  To  run  the  code, 
type  RUN.  The  code  will  iriterogato  for  pertinent  phy¬ 
sical  parameters.  A  sample  session,  run  on  the 
University  of  Tennessee  Space  Institute  VAX-.'i  1 ,  fol¬ 
lows  the  listing;  the  file  containing  the4  listing  was 
named  DEB 2 . BAS . 


PROGRAM  LISTING: 


:i  o  r 
'?  0  ! 
DO  f 
4  0 
HO  ! 

<;  o  r 

70  I 
no  i 
o  ;>  i 

:l  0  o 
1. 1.  o 

:i 

l 'AO 
,!  0  0 
I  O  ') 

1  00 
1  7  0 
mo 
:l  eo 
7 0  0 

:  o 
- ' '  •  ■:> 
N.'/O 


U  T 
I.)  T 


FH 
FM 
F  M 
!:"M 
FH 
KM 
I  H 
KM 

RTi'lf 
IMP 
I  HI  ' 
INPUT 
I  ’O  UT 
INPUT 

pp  ru  r 

REM  — 
REM- 
WF.II-  ~ 
REM" 
RE  M  • 
lif'il 
PPM  • 
REM- 


HE I’iRIG  IIAZAED  area  for  rocket  backblas  f 
This  r'f'odf  sin  ea  1  cm  J  c.  tee  the'  d.i.  mens  ion:;  of  &  circular  ru 
to  the  rear  of  &  rocket  .in  which  there  Is  ;?  debris  hazard 
1 1 'put  i  c  :i  i  •'  L  h  t*  c  e  r  ( ■?  c  t  i  on  s  {  o  n  e  s » ■  p  c  i  f  y  in  si  t  'n  ©  r  o  c  k  ©  t »  c  in. 
«rcci  (’Mind  object  or  >>c  r&on  In  bo  protected*  end  one  do  so  r 
the  debris.  The  prod  rain  is  designed  to  he  run  in  MKi*!  uni  I 
out  can  be  elnui  'oij  to  ane  other  system  of  con* is tent  uni  Is 
be  ehoruSirwJ  the  sir  ems  ton  t  «  in  eta  to  inert  to  330  end  340. 

*  ROCKET  BATA  * 

'  »  Tint 
"  >  Brio 

*  *  Pc 

*  >  K 

*  f  N  n  7.  a 1 1  ;| 

*  U  K  T  M  P  E  N  V.  T  R  A  T I U 1 1  P  A  f<  A  H  ETER  ' 

U  i  ,»  i’i  o  usurp  of  one  niy  i‘er  volume  t-ooui  red  to  fracture 
f.<r  lot  tint  t"  I-  i,  3  )  .  Tri  MK8  units  the  following  values 
veil  b<*  uio.’d  (  son  J.H. Lewis  <?t  a  1 1  'An  Emr-i.  r  ical/Mathcmatic 
Model  to  Estimate  the  Probability  of  HI*,  i  n  Ponot  ration  be 
V  a  r  :i  n  1 1  P  ro.it*  c  t  i  1  <  *  s  7  •  A I C  5 1.  T  R  -  7  G  0  0  4  *  A  j*  r  7  U  > 

boro  sKiii"  . 0  ■■  1  ‘J  0  0  0  0  k  >1  /  e  c  *  « 

2  I  a  y  o  r  i. 1  n  t.  f  o  r  in  -*  -  ly  -  A  4  0  0 00  f  >:i  /  •»  <?  c T  <:>  <„  e 
A  miner  un  i  f  m rm-~0  •7.,.0000  kd/suc$st»r 


t  O  '(' 


,i  c>  i,  n 

ft  y 


Nozzle  throat  diem 
H  o  z  z  I  e  c  j  t  d  i  >.•-  in 
Ch an. ho r  r rest: * i ro 
ih  eci  fic  t  ratio 
Nozzle  half  a  relict 


o.l. 


Best  Available  Copy 


'MO  ['JPUT  *  Q :  ‘  1 0 

■T  O  PRI  NT  '  DFBF*  \B  ji.VTh  Co nlei'  0  to  exit)* 

260  INPUT  *  Pi  cJDrljd  jpp3  1  *  >  fir 

270  REh - E;;  i  i  >  r  0 '  S'  J  ill  if  0  -iss  i  snetJ  to  HP. 

200  IF  AP- -O  Ti-S’H  GOTO  700 
2  7  T  N  i ;  r  “  Yi  3 «  15  1  *  f  M  P 

3 0 0  1 » P Uf  *  Bre 4  c o *-• f f  i e i e n it  ‘  ? C d 

3io  i\EM~* Phooir-densj.  ly  of  sir  in  kilosrpsns  r-cr  cubic  meter 
320  RJTh — P?j  »*=f'r?sfru,.,e  of  si5*  in  newtons  per  souf  re  nteter 
330  Rhosi r • t .2929 
:.‘J0  Pei  r;'  100000 

320  rxi-:— Lmult  is  lh&>  number  of  Jet  primary  y-ive]  ensths  ov;-r  which 
360  K L H  —  sistnif  itnjnt  particle  :<c  cole rsiion  occurs  *  Experiments 
370  REM---  c  v;  surfdesi  in  si  o  vc,  lue  other  thsn  1  is  more*  or-rr o?r  isia 


i  i ,  i  4  - 1  ' 

i  —  I  •  A  \.  S 

REM--IT1  *  !,2  ?  K3  Are  internal  constant 


p  r  o r  n  u. ;» i  n  a  c  on"  e  n  i  e  n  c  e  • 


•  400 

h*‘  1  -  }•’  J.  1  \  /  4  iv*  ^  1  *, 

P\i  \  I'll  A  v  .  i|\  I  / 

410 

!<  2  ~  2  /  f  K 1 

420 

K  “  2  (  !•■  "  1  )  /  (  J\  -i  1 

430 

K*  ii  i‘i  *■  **  h  r  ^  t  i  Q  ?■  s 

i  440 

{-  r  s  t  i  c  -  <  0  n /Bn  l 

4  SO 

REM- -The  nsxt  f 

:  ■  4  6  0 

REM---  Hsch  is  3 

!  t  4'0 

1 1 1 0  n  “  a 

Li  4&0 

ih-  ii  l=M;ich 

O  ;  \ 

i  i .■  ii  ~  S  Q  P  <  it  3  %  ( A  r 

il  5?b 

IF  ABS<Hoch-:{ex 

Li  53  0 

'130 

KCM*-*‘C '1*  u3  si-5 

n  ssq 

i:’  n ;;  i  fc  P  c (IT,  7  '}■ 

U 

R7M--Csl.  nl  r<ir 

^  1-1“  ^ 
O  *0 

U  ?  v  <*  .1  ~  1 6  9  0  B ;  1 1- M 

060 

REH-'-Cs  1  :.ulf.-ie 

570 

L-w3veliL.Mjii 

4  *  r..qo 

REM — Inl--rri3i  v 

5  7  0 

Porml  ~  •  5*i'-p-/  i  Co 

4  r  • 

|  l.*  V 

Ptrr.  ?  -i_Utt  <  2  i-LTC 

II  /  -»  A 

;  |  ^.11 

3%  L.  i  1 - L  i  C*  Ui-il*? 

320 

3  L  ■'  r  *  s  0  f  ,  '■  P  3  r  t  >!' 

n  6  3  *»• 

1'  O  i  Hj.  n  .1  D 

H  : 0 

PCM . R 0  m  d i n  0  2 1 

REM- -Cslculclo 

—  ^60 

M.-ir'-S0R(K2:-:.<  (P 

1  !  /  /> 

1  i  *  •  ’ 

REii--Cslc.il  pic 

1  '■  6  P-0 

f 1  h  I* r*»  /•  ■  •  H »-  '  j  1  3  **  i 

690 

Ar.Mlo3--M'0R<Ki'-x 

1  i  700 

t  1  »-  “ ‘Ur  3  A  M  •  J  A  c  !•  ? 

!  i  710 

Psrn-3  M::  ir '2-1 

7  3  0 

OnC  r-7  r.'  jR  ( !%  J  )t 

n  730 

!\  £  H  --  ■  C  3 1  ■'  u  1 :  t  <y 

1  1  740 

1  ! 

B  i  s  r--  c* r.ii-57  s  3 f  f  A 

•■SO 

PkIMT  *  uur.? 

„  730 

PRINT  '  * 

1  | 

fimo  ;  00 

biPrtdoi  1 


"  ■  S  U  K  •  F  3  r  r.  6 ,  K I  >  >  -  ft !  ,N  ‘  F.  &  R  <  P  3  ;•  n.  3 )  > 
sndil -Hewer  Cumins!  snalc  for  Hair* 

t'(  C  ;•  G  K  \  l  '  {  IT,  y  ,  l\  j  )  /'"Si  I  |<E  (  i'UR  \  P*J  I’  :  ■  1  I  ) 

s:  e r  cion  hr}  f  civile  f  Risi-ansif  in  degrees, 

1  v  2  -  »*»  n  U 1  c  1  )  -=  N  o  7  3  n  s 

ion  up  io  '  » Bit? end  -  *  decrees* 
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Figure  E.l  -  Example  session  of  BASIC  debris  hazard  code 


i.*  iibii 

iX-1  j 
R  c?  b  d  y 

Hi!  Tt  nppn  BAG 


5-KEC  -1980  17: 18 


i 

I 

- 


Mi 

r: 


di  3.-:!  J 

?  *  0424 

V*  -t  3  !it  « 

7  ,0616 

3  *•>  u  i  : 

?  6 B000000 

ratio: 

?  1,16 

pt-.:0  <?  { 

?  11 

PARAMETER 

q  : 

?  340000 

l er  0  to 

exit  > 

b  r  e  a  5 

?  ,0001 

Hass  { 

?  ,0007 

cisnt t 

?  ,3 

standoff 

di&tsr.cc  = 

to 

77.3552 

ter  0  to 

exit) 

t-  S  i 

?  0 

33,6709 

decrees 
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B AC KB LAST  DEBRIS  HAZARD  CODE  FOR 


HP-lllC  PROGRAMMABLE  CALCULATOR 


PROGRAM  MOTES : 


I  a 


'i  I 


This  code  requires  63  registers  of  program  memory 
and  15  registers  of  storage  memory.  Running  the  code  puts 
the  calculator  into  degrees  mode.,  fixes  the  display  at  zero 
digits,  and  sets  flag  01. 

Lines  ,001  through  03 are  for  input;  line  025 
is  an  entry  point  for  debris  specification.  Lines  037 
through  068  solve  equation  (2.22)  for  Mach  number  as  a 
function  of  area  ratio;  MF  remains  in  register  015  after 
the  program  runs.  Lines  070  through  0-95  calculate  the  jet 
wavelength  and  store  it  in  register  03;  1  remains  in  that 
register  after  the  program  runs.  Lines  097  through  12'4 
calculate  the  standoff  distance.  Lines  127  through  1A6 
calculate  the  dispersion  angle  0  Lines  1^7  through  15^ 

display  standoff  distance  and  maximum  dispersion  angle. 

Line  156  loops  to  additional  debris  input.  Lines  157 
through  17^  are  a  subroutine  for  calculating  the  Pranatl- 
Meyer  turning  angle,  equation  (2.31),  given  Mach  number. 

The  code  is  set  up  to  handle  multiple  debris 
specifications  consecutively.  After  the  first  specification. 


flag  01  causes  substantial  portions  of  the  code  to  be  skipped. 
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PROGRAM  USAGE 

The  listing  below  is  designed  for  SI  units  arid 
for  standard  atmosphere.  Line  005  has  ambient  air1  pressure 
of  100,000  N/m2  and  line  120  has  ambient  air  density  of 
1.25  kg/m  .  The  code  can  be  set  to  ether  units  or  to  other 
ambient  conditions  by  changing  the  numerical  constants  on 
those  two  lines.  Data  entry  must  be  in  consistent  units. 

•  To  run  the  code,  load  and  press  R/S.  An  example 
session  is  given  in  figure  (E.2). 


PROGRAM  LISTING: 


001 

LBL  T DEBRIS 

016 

TCH AMBER  P 

002 

CP  0 1 

017 

PROMPT 

003 

DEG 

018 

STO 

06 

004 

FIX  0 

019 

tGAMMA 

005 

100000 

020 

PROMPT 

006 

STO  01 

021 

OrpA 
oi  J 

07 

007 

tTHR0AT  diam 

022 

T  PEN 

ETRATION 

00  8 

PROMPT 

023 

PROM 

FT 

00S 

STO  03 

024 

STO 

08 

010 

tEXIT  DIAM 

025 

LBL 

01 

011 

PROMPT 

026 

7  DEBRIS  AREA 

012 

STO  0-4 

027 

013 

rN0ZZLE  l 

028 

STO 

09 

014 

PROMPT 

029 

T DEBRIS  MASS 

015 

dll'  JS J 

030 

PPAM 

i  iiv;  s 

prp 

085 

RCL  01 

086 

/ 

087 

SQRT 

088 

RCL  04 

08s 

& 

090 

RCL  15 

091 

SQRT 

092 

s 

093 

•  69 

094 

s 

095 

STO  03 

096 

LBL  04 

097 

RCL  03 

098 

STO  04 

099 

BCL  11 

100 

$ 

101 

RCL  06 

102 

s 

103 

104 

¥. 

105 

p-pr 

*  -  v/JL?  pg 

106 

/ 

/ 

j 

T  s'* 

j  ;]J 

108 

RCL  07 

109 

SQRT 

110 

T  SI 
!-»!« 

112 

* 

113 

•!- 

114 

RCL  10 

115 

a 

116 

RCL  11 

117 

/ 

118 

RCL  09 

119 

/ 

120 

1.29 

121 

/ 

122 

2 

123 

/ 

124 

ST  +  04 

125 

PS?  01 

126 

GTO  06 

127 

SP  01 

128 

RCL  06 

129 

RCL  01 

130 

/ 

131 

1 

132 

RCL  07 

133 

1/x 

134 

- 

13? 

vx 

lju 

1 

137 

I  II 


Figure  (£.2).  Example  session  of  HP-^1C  Debris  Hazard  Code. 


USER  ACTION 


DISPLAY 


REMARKS 


R/S 

THROAT  DIAM 

Enter  nozzle  -throat 
diameter,  d* 

.0*12*1 

.  0^12^1 

d5  in  meters 

R/S 

EXIT  DIAM 

Enter  nozzle  exit 

diameter,  d«  . 

.  06l6 

.0616 

d_  in  meters 

H. 

R/S 

NOZZLE 

Enter  nozzle  diver¬ 
gence  angleOv, 

- 1 

11 

11 

G,.  in  degrees 

R/S 

CHAMBER  P 

Enter-  chamber  pres¬ 
sure  P 

O 

680C00G0 

68,000,000 

PQ  in  N/;n“ 

R/S 

GAMMA 

Enter  specific  heat 
ratio  k  of  exhaust 
products 

1.16 

1.16 

R/S 

PENETRATION 

Enter  penetration 
parameter  Q 

P 

0  U 

JTvUvU 

3*10,000 

Qp  in  kg/s2 

R/S 

DEBRIS  AREA 

Enter  particle’s 
projected  area  A, 

6 

.0001 

.0001 

A  in  m 

R/S 

O  TO  M  A  Q  C 

Enter  particle's 
mass  m 

.0007 

.0007 

m  in  kg 

R/S 

DRAG  COEPF 

Enter  particle's 
drag  coefficient  C- 

APPENpIX  E 


/  t 


SOLUTION  OP  EQUATION  (2.22) 


Equation  (2.22)  gives  area  ratio  as  a  function 
of  Mach  number: 


Oc+1) 

2 f k-l) 


(2.22) 


The  model  requires  Mach  number  as  a  function  of  area 
ratio.  The  inversion  can  be  accomplished  by  iteration. 
Reformulate  equation  (2.22)  to  read: 


i 


\ 

/ 


Assume  a  trial  Mach  number  of  unity  and  evaluate  the  right 
hand  side  of  equation  (E.l).  Take  this  value  as  an  updated 
trial  Mach  number  and  repeat  the  procedure.  Satisfactory 
convergence  will  generally  occur  in  four  or  five  iterations. 
See  also  figure  (3-3). 


